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SYNOPSIS 


In "the present work two dimensional harmonic wave 
propagation in fiber reinforced composite materials has been 
attempted using Finite Element Method (FEM). Concepts of periodic 
unit cell and Floquet’s conditions has Vjeen used in the analysis. 

Wave propagation analysis is known to give a good indication 
of static and dynamic behavior of the medium of propagation. 
Nature of wave propagation in isotropic materials has been 
extensively studied by many researches and exact results for the 
same has been obtained. Eecently there is increase in the use of 
composite materials in critical and non critical components like 
aerospace structures. Understanding the behavior of these 
materials under different conditions such as impact and fatigue 
loads becomes important. Non destructive testing techniques like 
ultrasonic and acoustic emission techniques which are used to 
study the material characteristics are based on the principles of 
wave propagation . As conventional techniques are not adequate, 
study of wave propagation in these materials becomes necessary. 

Attempt has been made to study the wave propagation in 
composite materials to obtain ten lowest eigen frequencies and 
corresponding eigenvectors as a function of wave length X . 

Two different composite materials ( E. /E = 3.275 and 

1 Ttl 

20. Dare studied and their behavior is compared to that an 
isotropic material. Frequency spectra of these materials for the 
propagation in the plane of fibers and in the thickness plane have 
been obtained. The wave patterns corresponding to optical and 
acoustical branches of the frequency spectrum are investigated. 

Results reveal that composite materials are dispersive in 
nature as compared to isotropic mamterials. The amount of 
dispersion depends on the wavelengbh/f requency , direction and the 
mode (acoustical and optical) of propagation. It has been observed 
that dispersion increases with the ratio of stifnesses of 
individual constituents and in the propagation perpendiculer to 
the fibers.lt has been found out that FEM can be effecitively used 
to study wav© propagation in composite materials. 



CHAPTER 1 


INTRODUCTION 

1.1 DIN IT IR<DtDq,)CHO (PIN : 

Static and dynamic analysis of composite materials is gaining 
much attention these days as these materials axe finding wide 
range of uses in many critical and non critical structures of 
aerospace structural components. These applications impose severe 
restrictions on the behavior of composite materials. Especially 
under conditions of high strain rate and impulse loading 
(shockwave during impact) their behavior is to be understood in 
sufficient detail. Non destructive testing (NDT) techniques like 
Ultrasonics, Acoustic emi.ssion and Acousto Ultrasonics are being 
widely used as experimental tools to .study the structure and 
hence, • behavior of materials both qualitatively and 
quantitatively. Also since there is a relationship between 
mechanical properties and wave propagation, study of wave 
propagation gives another method of determining the mechanical 
properties apart from conventional methods, like static and 
vibrational propertie.s. 

1.2. tHAS/E F[R<DFA€AirD<D[N ; 

Harmonic wave propagation analysis deals with studying the 
mechanics of wave propagation like stress wave propagation in 
materials. This analysis to determine properties such as E, 
has got certain basic advantages over other methods of analysis, 
Many complexities that are involved in the solution of governing 
equation of motion can be avoided if the wave solution technique 
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Is applied to such problems. Though modal superposition analysis 
can be used to obtain free/foroed response of continuua they are 
not often applicable for cases with complex geometries, plates or 
systems of infinite or semi infinite exbent. Sometimes it may be 
difficult to obtain mode shapes or frequencies for medium v?ith non 
unifoimi properties and heterogeneities. A medium of infinite 
extent has a continuous band of frequencies and hence mode 
superposition etc., may not be directly applicable. Some of the 
above difficulties can be circumvented by applying wave 
propagation technique, 

1.3. CC>lh[POSDTE mATTERDAn-S ; 

A composite material consists of two or more matez'ials mixed 
in suitable proportion so as to obtain desired qualities. Usually 
the proportion is determined based on the optimal xitilisation of 
properties of the constituents. There will be generally two kinds 
of constituent materials, viz., a reinforcing or the major load 
bearing material and a bonding matrix ma-J>erial. Composites are 
generally classified into particulate reinforced composities and 
fiber reinforced composites depending upon the geometry of the 
reinforcing material. The high strength and stiffness, high 
directionality in the elastic properties (anisotropy ) and good 
formability are some of the attractive features of composite 
materials. 

The complex nature of composite material presents 
many difficulties in analysing composite materials and to 
tinderstand mechanical behavior due to the heterogeneity at the 
micro level ( when the area of Interest is of the order of fiber 



3 


dimensions ) because of differenb material properties , 
interfacial effects between plies, and viscoelastic behavior of 
the composite. The general anisotropy and coupling of in plane and 
out of plane stresses / strains introduce lot of complexities in 
the analysis. 

The material has some inherent defects such as fiber breaks, 
matrix debonding, dela7ninatiori etc. . These cause adver. 5 e effects 
in the behavior of these materials under different types of loads. 
Plastic and viscoelastic nature of matrix material cause flow of 
the material beyond a certain value of stress. These factors have 
lead to Increasing use of wave propagation analysis as a technique 
to study the structure and behavior of composites under different 
loading conditions. Standard methods and solutions for isotropic 
materials cannot be easily extended to the case of composites. 
Various theories have been proposed to model the elasto dynamic 
behavior of these materials. 

1.4. ITDlHD'iriE EILE{HE1N¥ [HE¥IH<P 0 > ( F E tH ) : 

It is an approximate analytical technique for solving 
differential equations. Using variational or weak formulation of 
the differential equation and then discretising the domain 
suitably a set of algebraic equations are obtainedDO. These 
equations which are easier than the original differential equation 
can be conveniently solved using a digital computer. Finite 
element method has several advantages over the exact solution 
techniques and other approximate methods. This technique can 
easily handle complex shaped domains , complicated boundary 
conditions, discontinuities, nonhomogenieties , nonlinearities in 
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material etc.. Advent of high speed digital computers has 
facilitated rapid increase in the range of applications of this 
technique in solving realistic problems. 

Fi'om the above discussion it can be inferred that the wave 
propagation analysis using FEM can be a good technique to 
characterise the wave propagation composite materials. 

l.S. OFElRAFqJIRIE : 

Literature available in the field of finite element analysis 
of wave propagation in composites is limited in comparison to 
that available in vibration and modal superposition studies. Some 
work is done in wave propagation analysis using finite difference 
schemes and special matrix methods [2,. 3]. Kohn et al, [5] used 
Rayliegh Ritz procedure to solve variational form of wave equation 
formed over a single periodic cell using Floquet theory and 
determined dispersion properties and elastic properties of a few 
composites using plane wave approximation for two and three 
dimensions. Natural frequencies and mode shapes obtained are 
comparable to the available exact results. 

Dong and Nelson [6] have studied natural vibrations and wave 
propagation in laminated composite plates using FEM. They have 
obtained ten lowest natural frequencies and associated modal 
displacements, along with stress distributions. 

Yang and Lee [7] have used Floquet theory and finite 
difference scheme for modal analysis. Lowest five modes are 
accurately and efficiently computed. 

Nelson and Navi [4,8] have used FEM to study harmonic wave 
propagation in composite materials and other orthotropic materials 
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under plane strain condition. Flocjuet tVieory and periodic cell 
approach has been folio-wed. They have shown that this model can 
represent accurately wavelengths upto half the lattice dimensions. 
Frequency spectra for various orthotropic and isotropic materials 
are obtained. The eigenvalue system has been solved for various 
values of linear wave number C . 

Golub et al. [9] have solved an eigen value problem obtained 
using FEM procedure. They considered steady state 2-D wave 
propagation in problem and used Floquet theory to obtain a 
Hermitian system of equatio'ns which are then converted to real 
system of equation using the sparse nature of the matrices. 
Dispersion relations are obtained from the analysis. 

Minagawa et al. [10,113 have used finite elements to study 
harmonic wave propagation in fiber reinforced and layered 
composite plates. Unit cell approach is followed here based on 
Floquet theory. Wave patterns corresponding to acoustical and 
optical branches (in the Brillouin sense) of frequency spectra are 
obtained for different wave numbers. It has been observed that for 
wave lengths greater than the periodic dimensions of the medium 
and for low frequencies, the medium behaves like a homogeneous 
one. 

Datta et al. [12] have used a stiffness method to study wave 
propagation in oomposi-te plates. They have found that frequency 
spectrum for ^ laminate differs considerably from that of a 
(0/90/0) cross ply laminate, 

Osaki and K^ropara [2] have made an attempt to study defects 
in unidirectional fiber reinforced composites by elastic wave 
propagation. They have used a finite difference scheme for 
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analytical modeling and ul^trasonic "tecbnique fox experiments . By 
defining a new parameter obtained using FFT analysis of 
experimental data they have evaluated the defects in the form of 
fiber breaks . Vector charts have been used to illustrate the effect 
of fiber breaks (notched defects) on sine wave propagation based 
on dynamic simulation. 

Mai and Barcohen[13] have, in their study of fiber composites 
used a matrix method for analytical simulation of leaky lamb wave 
(LLW) experiments . They have calculated phase velocities of Lamb 
waves of laminates. Ref lection and water loading effects have also 
been considered . For simple cases they have found good agreements 
whereas for realistic cases discrepancies have been noticed. 

Toshiuki Oshima et al.[14] studied the effect of the 
orientation of point masses on stress waves. Dynamic response of a 
rectangular cross sectioned beam with regularly arr-anged point 
masses under impact loading is investigated. Results are compared 
for the cases with discrete and continuum models. 

Masanori Koshiba et al, [153 have used FEM and an analytical 
method to study Lamb wave scattering in an elastic plate wave 
guide. Effects of internal and surface cracks are investigated. 
They have found out the effect of wedge angle (of the crack) on 
the reflection coefficient and associated resonance phenomena. 

Considerable work was also done to characterise composite 
materials using NDT techniques. A brief report of the same is 
provided here. 

Rose [16] has underlined the principles of ultrasonic wave 
propagation in composite material inspection. The effects of wave 
velocity, dispersion, ref lection coefficients and attenuation 
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charactei'istics of the material on the wave profile are studied. 

Rose et al. [17] attempted to analyse actual wave profile in 
an unidirectional graphite epoxy plate u-sing through transmission 
technique. Wave velocities are measured for different fiber 
orientations , From the wave velocities^ Young’s modulii in 
longitudinal and transverse direction are computed. Other 
engineering constants were also computed. They reported good 
agreement between the values obtained from experiments and those 
predicted by the theory. 

Tauchert and Quzelsu [18] have used ultrasonic pulse 
technique to study dispersive I’elations of plane harmonic wave in 
boron epoxy composites. Acoustical velocities in longitudinal and 
transverse directions are measured for different frequency pulses. 
The relationship between the group velocities and frequency 
revealed the dispersion characteristics of the material .. The 
static modulii for the composites were also computed. 

Munson and Schuler [19] proposed a mechanical theory to study 
the wave propagation in laminated media. Both longitudinal and 
transverse propagation are considered, linear elastic models and 
nonlinear hydrodynamic models are proposed for different materials 
and pr©s.5ure levels. They have found out that propagation normal 
to the laminates is identical to propagation along the laminates. 

Hemann and Baaklini [20] have used ultrasonic technique to 
show that stress wave propagation is dependent on stress 
levels . Attenuation of stress waves was found to be related to 
tensile stress levels in the composite. The stress levels were 
kept within the elastic limit and no macro level damage was 
observed. Wave speeds in the graphite ©poxy material were found to 
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vary between 2700 m/s and 6800 m/s for pulses of frequencies from 
15 kHz to 1.9 MHz. It was ob.5erved that wave speeds increase 
slightly at all tested frequencies with increase in stress levels. 

Recent developments in digital signal processing and 
electronic hardware engineering has facilitated better extraction 
and analysis of experimental data and thus has enabled application 
of the same for defect identification and quantification. 

Teti and Caprino [21] have used digitized ultrasonic data for 
analysing defects like porosity delamination etc., in thick glass 
fiber composites. Complete waveform digitization and analysis has 
enabled them to localize and to identify the defects. 

Dayal and Kinra [22] developed a new technique for ultrasonic 
evaluation of thin specimens. Use has been made of FFT methods 
with conventional ultrasonics for measuring phase velocity, group 
velocity and attenuation of thin (sub millimeter and .sub 
wavelength thick) specimens. The technique has been found to be 
suitable for dispersive materials also. Effect of number of 
transverse cracks on the attenuation has been studied in graphite 
epoxy laminates. Attenuation seems to be increasing monotonically 
and significantly with increasing nun^ber of cracks. Also phase 
velocity seemed to remain constant with increase in damage. 

Ultrasonic technique has been used to study the fatigue 
behavior of composites [24,25]. Variation in time and frequency 
domain parameters of the output pulse have • been studied to 
associate them to specific changes in the structure and 
characteristics of the material. So far one to one correspondence 
between the two has not been established conclusively. 

From the literature survey it is evident that much of 
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analyiiical work has been concentrated on studying dispersion 
relations and modal frequencies and displacements in composites. 
Main difficulties that arise while analysing these materials are 
their dispersive nature, heterogeneities at the micro 
level, viscoelastic effects of matrix etc. Some Work has been done 
to overcome these difficulties using alternate models such as 
control volume analysis C26] etc. Application of effective 
stiffness and effective modulus theories to study dynamic behavior 
of compo .sites are reported in references 27 and 28. Both discrete 
and continuum models are considered . Tewary [27] has extended the 
same to study the effects of defective fibers on wave propagation. 

1.6. PIRESEfNir IHOIRIK : 

In the present work attempt has been to study the harmonic 
wave propagation in composite materials. Emphasis is on 
determination of dispersion relation in the plane of fibers and in 
the thickness direction. Analysis is done by two methods. 

(a) By modelling the composite material in a such manner that 
fibers are uniformly spaced 

Cb) By modelling the composite medium as a uniform homogeneous 
medium with appropirate average properties. 

Concepts of Floquet’s theory and finite element method ar-e 
used in the analysis. Dispersion relations are determined for 
various material properties. 



CHAPTER 2 


BASICS OF WAVE PROPAGATION AND ITS APPLICATION IN COMPOSITE 

MATERIALS 

In this chapter a brief overview of the underlying principles 
of wave propagation and basic analysis procedure of composite 
material is given . 

PlR<DlPAC&'[rD«P!N : 

Wave or disturbance is a process or an influence like some 
form of energy which moves through the material with a finite 
velocity, without causing any bulk transfer of the material as a 
whole. 

Waves may be either standing or travelling waves. In 
vibration problems one encounters standing waves, such as in the 
case of vibrations of strings. In case of a shock or impact 
loading the travelling waves are generated . Depending on the form 
of the solution desired for the wave equation one gets a standing 
wave or a travelling wave solution. Separation of variables 
technique yields a standing wave solution whereas the method of 
characteristics gives a progressive wave solution. 

The governing equation of wave motion is of the form 
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where is some disturbance, c, is a physical constant and Xj. the 
coordinates, in x,y,z directions. Choosing v'(x,t) = X(x) 

T{t) leads to a standing wave solution and V' = V'(x-ct), gives a 

travelling wave solution. 

2.1.2 IH&KDiOMOC ; 

If the solution is of the form: 

( 2 . 2 ) 

where X is the wave length, c the velocity of propagation and 
is the wave number giving the number of waves per cycle , a 
harmonic wave solution is obtained. The relation between wave 
number k , circular frequency « and velocity of propagation a is 

w = kc (2.3) 


V = A sin 


[ 


2rr 

~ 


(x + ct) 


2.1.3 tPmSE AfNiD <ElR<D<UtP yEtt-cpCfllTDES : 

The velocity with which the points of same phase move in a 
medium is called as phase velocity. This is same as c in equation 

2.3 . 

The velocity with which the modulations or beats propagate in 

the material is termed as group velocity. It is denoted by c and 

g 

is given by the relation 

do> 

c ~ (2,4) 

® dk , 

nodes are patterns of motion having a property that at any 
point in. it; the body moves perfectly sinusoidally and all points 
move with same frequency. 
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Dispersion is the frequency dependent phenomenon of a medium 
whei'ein the wave front undergoes a change in shape as it 
propagates in the material > Flexural waves and waves in 
anisotropic material are certain examples of dispersive waves. 

2.1.4 irW[ES: <DF [HAVES : 

There are different types of waves called longitudinal waves, 
transverse waves .plate waves, surface waves etc. In longitudinal 
waves the direction of particle motion is same as the propagation 
direction. If the direction of particle motion is perpendicular to 
the direction of propagation then the wave is called as shear or 
transverse wave. Sound waves travelling in air are examples of 
longitudinal waves and electromagnetic waves are examples of 
transverse waves. 

Waves are also classified a.s bulk waves , sjjirf ace waves etc.. 
Longitudinal and shear waves travelling in an infinite medium are 
examples of bulk waves. Here the surface effects are not felt by 
the waves . 

Surface waves are waves propagating near and on the surface 
of the medium. A plane surface wave is called a Bayliegh wave in 
which the particle motion will be elliptical. Waves caused by the 
movement of earth’s crust and travelling on the surface are 
examples of Rayleigh waves. Surface waves travelling in thin 
plates and shells are termed Lamb waves or plate waves, 

2.1.5 nCxaHAEF TT[HIE«DIR^ : 

Floquet theory concerns with the solutions of differential 
equations with periodic coefficients. Consider a differential 
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equation of the form , 

a^(x) y (x) ^aJ^(x) y (x) + ag (x) y(x) = 0 (2.5) 

in which the coefficients complex valued, piecewise 

continuous and periodic, all with tlie same period 'd': 
i.e., a^(x+d) = aj.(x) , { 0 < r < 2) and ’d’ is a non zero real 

constant. It can be observed that if V'(x) is a non trivial 
solution of equation 2.5 then v'(x+d) is also a solution. Floquet 
theory . provide.^ a proof for the property that 

V^(x+d) = p v'(x) (2.6) 

where p is a non zero constant. 

This result finds application in solution of the wave 
equation and it simplifies the solution appreciably, because it 
facilitates application of concept of vmit cells and periodic 
boundary conditions to the problem. 

2.1.6. <UNDir CED-B, A(P(F(R<DAaH : 

As mentioned above Floquet theory simplifies the analysis of 
the wave problem in a periodic medium to the analysis of a small 
a portion of the spatial domain, called a unit cell which 
satisfies the following conditions. 

1. The medium under consideration can be discretised into 
identical cells representing one period of the medium, 'd’ . 

2 All the characteristics of the medium should be represented by 
the cell* 

Fig. 2.1 represents a one dimensional medium discretised into 
unit cells. If it is assirmed to contain N number of particles in 
each cell, the particle N+p is said to be related to particle N+q 
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Pig. 2,1 Unit cell in a periodic one dimensional medium 



Fig, 2.2 Unit cell an isotropic 
3-D infinite medium 
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Fig, 2,3 Unit cell in an^infinite 
3-D infinite iieaium 
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in the same way as particle p is related to particle q. It can be 
seen that if N approaches infinity then the discrete medium 
becomes continuous . This concept can be extended to two and three 
dimensions also. For a general three dimensional medium the 
representation of unit cell is shown in the figure 2.2. 

In a homogeneous material there is no physical significance 
for an unit cell. For homogeneous crystalline materials the 
lattice can be lil^ened to an unit cell. In amorphous materials it 
is difficult to identify a basic representative and repetitive 
unit. 

A uniform composite material can be idealised as the fibers 
being arranged in regular array and embedded in a matrix material. 
Here the unit cell can be seen to be a region shown as dotted 
rectangle (Fig. 2. 3). In some cases the unit cell may be hexagonal. 
It is evident from the figure that the cell pos.sesses all the 
engineering properties of the medium . In the present study 
rectangular type of unit cells are assumed to discretise the 
medium. 

2.2. {HAVE ClHAtRACTrElRDSTrDCS fl(N COiHIP€>SD¥E lhft7EIRD/&0-S : 

Wave characteristics in homogeneous materials have been 
studied extensively [29,30,31]. But only limited understanding of 
the same is available in the case of composite materials. Material 
non homogeneity, layered structure , anisotropy etc,, that 
differentiate the composites from isotropic materials are 
reflected in their wave propagation characteristics. Figures 

2. 4, 2. 5 [16] show the difference in wave propagation 

\ 

characteristics of the two types of materials. 



ir 




a 


b 


Fig, 2,4 Spherical wave velocity profile for a 
homogeneous isotropic material. 


' Fig, 2,4b A typical wave velocity profile for a composite 
material , 



Fig, 2,5a Typical A-mode display for a homogeneous 
isotropic material. 

Fig, 2,5b Typical A-mode display for a composite 
material , 
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Composites materials are dispersive in nature due to their 
inherent anisotropy and non homogeneity. Thus the wave profile is 
constantly changing and it is difficult to separate different 
frequency components of the wave. This also implies that it is 
difficult to define phase velocity for these materials. The 
velocity that can be measured is the group velocity. 

The acoustic imi>edance . represented by pc is difficult to 
obtain in the case of composite materials as c is a function of 
frequency. Since this property is important in assessing the 
propagating characteristics it is often necessary to define an 
average value of the acoustic impedence for specific cases in 
order to study the characteristics. 

When a wave passes througVi an interface of two media some 
portion of it get reflected and some get transmitted (refraction). 
The refraction formula Is also known as Snell’s law^ written 
mathematically as 

Cg Sin0^= Cj^Sin^g (2.7) 

where ©j^and ©g are shown in Fig. 2.6. This cannot be easily 
applied to composite materials as velocity is dependent on the 
direction and the associated mode conversion . 

Attenuation levels are much higher, in composites than in 
Isotropic materials. This effect is more pronounced at higher 
frequencies as the associated wavelength approaches the size of 
the fiber , Local heterogeneity is also a factor that causes 
attenuation due to scattering of waves ebc. . 
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5 h 


Fig. 2.6a, Refraction of waves at an interface between 
two isotropic media. 

Fig. 2,6b, Refraction and mode conversion of waves at 
an interface of a composite medim 



CHAPTER 3 
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PROBLEM DEFINITION AND FINITE ELEMENT ANALYSIS 

In this chapter the present problem is defined and necessary 
assumptions made are described. Finite element formulation of the 
problem and its analysis procedure are highlighted. 

3.1. FtROELElM DEFDlNnFDOEN : 

In the present work attempt has been made to study wave 
propagatio7i in composite materials. Attention has been focused on 
two dimensional harmonic wave propagation in fiber reinforced 
composite materials. 

In two dimensional wave propagation two cases are of 
interest vis. , the propagation in the longitudinal plane 
containing fibers and in the thickness plane as shown in figiare 

3.1. The waves propagating in these two planes can be either 
parallel or circular/elliptical. Bence the cases considered are 

a) Wave propagation in plane of the plat© (Fig. 3.1) with the 
particle motion confined to this plane alone. 

b) Wave propagation in thickness direction (Fig. 3,2) with the 
particle motion confined to this plane alone. 

Following simplifying assumptions are made in the analysis: 

1. The composite material is defect free and has fibers 
distributed uniformly. 

2. The continuum can be represented by unit cells as explained in 
section 2.1.6 

3. The material behavior is to be elastic in nature. 





4. The medium is of infinite extent in order to avoid 

boundary effects such as reflection and refraction. 

5. Case (a) is treated as plane stress case and case (b) as plane 
strain case. 

The periodic nature of wediium is illustrated in figure 2.3 
and typical unit cell. Because of the periodic structure of the 
medium Floquet’s theory can be applied giving a simple and 
feasible method of analysis. 

The governing equation of motion is of the form 


1 Fc. r 1 1 

2 5x. [Hjvi [ J J 



(3.1) 


where p is mass density , is elastic stiffness coefficient 

3) and u j^are generalised displacements (u,v,w 

components in x,y ,2 directions respectively) . The displacements are 
assumed to satisfy the harmonic relation , 

uf,y.t) = u^. 7 )e (3.2). 

Which results in an eigen value problem. This is solved for 
natural frequencies, co and associated mode shapes (eigen vectors) 
as a function of wave number C . 

3,2. FDD'IDFE ELElhlElNTr r<&tR(HaJa.AS>irDCHH : 

The wave equation 3.1 is solved by FEM u.sing the Galerkin 
method [1]. Assuming an approximate solution for u^^ in equation 
3.1 as Uj^and applying weighted residual principle, the residual 


can be written as 
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R C v>^ .V ) = 




P dA (3.3) 

dt ^ 


where w =w(x,y) are the weighting functions and R ,tVie residual 
over the domain D. Considering the first term in E and applying 
divergence theorem , it takes the form: 



The first term in expression 3.4 is the boundary term and as 
along the boundary are assumed to be zero, the residual R can be 
written as 






at' 




dA 


(3.5) 


Using equation 3.2 in the residual , R Eqn.3.5 takes the form : 

R = J I [ ^ ^ 

D 


3 . 3 . FD1N07E ELEWEiNI 01 )SClRIETrDZAirfl<l>lH : 

Figure 2.3 show.s the domain and a typical unit cell, the 
origin of which is located at n, + n.,d„ from the system or domain 
origin. An arbitrary cell is discretized into rectangular 8~noded 
isoparametric elements (Fig. 3.3) . The generalized displacements 
at every node are u and v in x and y directions respectively, 
satisfying continuity. 




We have 
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V = E N 

tv 

u = E u N 

LI 

V = E V. N (3.7) 


where N, are the nodal shape functions , in natural coordinates K 
and f) (figure 3.3), and i the node number. 

if ••=1 (1 + K ?,)(!+ h h.)(- 1 + ^- ?. + W. ) 

V. *4 V V It 

i = 1, 3,5,7 

if = I (1 - + h T). ) 

i = 2,6 

if - I (1 - n^'Xl + ? ?. ) 

i = 4,8 (3.8) 

The nodal displacements u^, are , for node; 


Ti Wnk + nk." wt ) 
u = U. e 2 ^ ' 

i ) 


V. = Y. e 
j i 


i(n^l5: + t^'t) 


(3.9) 


In equation 3.9 ^ 2 ^ x'epresent origin of the unit cell in 


term.s of wave vector. If U. and V. are expressed as follows, 


U. = Ul. + i U2. ; 

3 i J 

V. = VI. + i V2. (3.10) 

j i j 


the nodal degree of freedom vector takes the forni; 



L Ul. U2, VI, Y2 Jj 


(3.11) 
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3.4 ELEIHEIHTTAO. [HAITIRDCIES : 

The nodal shiffness and mass mat-rices are of "the form [32]; 

T 

[ K 3. . t J [B], [D], . [B], dA 

D 

C m ] = J p inf [N], dA 


D 


(3 


where [B], is the strain matrix and [D] is the coefficient matrix 
of the form'- 




0 


0 



0 


N. 


<■ .y 


0 

N. 


,y 


N. 

»■ ,y 
N. 

J .X 


N. 


>- >y 
N. 


J ,x 


[D],. = 


D. 


Sym 


D_ 


D 






D, 


s> 


(3.13) 


From equations 3.12 and 3.13 nodal stiffness and mass 
matrices are 4X4 and elemental mass and stiffness matrices are 
32X32 in siae. 

The integrals in equation 3.12 are evaluated using numerical 
integration . Using Gauss quadrature integration rule they can be 
expressed as : 

[K].. = 5: z W W [B3. [D]. [B] . |J| 

[M].. « E Z W W p N.N.|J| (3.14) 

ij p q p q V J ' ' 

where W and W are the weights of the quadrature, p and q are the 
p q 

number of points in ^ and n directions and jJ} is the jacoblan of 


. 12 ) 



transformation written as; 

|<^| = ( x.j, y,^ - x.y^ y,^ ) (3.15) 

It can be seen from Eqn.3.14, 2x2 integration is necessary 
for stiffness matrix calculation and 3x3 integration for mass 
terms. The summation is to be performed for 4 times for every 
stiffness term and nine times for every mass term. The special 
nature of [D] can be made use of in order bo simplify 
computation. Using an efficient algorithm [32] one can foi'm 

an intermediate matrix containing the summed up products of shape 
function dex'ivatives , of the form ; 


[A].. = 


2 S N. N. 

V ,X J,X 


2 2 N. N. 
^ )y 


2 2 N. N. 

^ J fV 

CA^) 


2 2 N. N. 

^ ,y j .y 

CA,) 


Wp W, |J| 


(3.16) 




It is enough if we compute [A]., by summation at four gauss points 

^ J 


and then compute [K]. . , thus saving apprecdable amount of 

^ J 

computation time. It is possible to simplify further the 
computation of [K].. using the special nature of [B]. . By forming 
another matrix [S], . written as follows, 

VJ 


CS],/-- 


CSll) 


(S12) 

(S22) 


(3.17) 


i j 


(S21) 
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[Kj^ . can be wribben as follows; 



Similarly for [M]^j given by Eqn. 3.14, one can form an 
inbermediabe matrix bhab stores the summed up product of shape 
functions, N and N and then express [M.].. as follows. 

^ si * y 



where L. . - £ Z R N. \J\ W V (3.19) 

tj t j ’ • p q 

p q 

Thus bo compute elemental stiffness and mass matrices CK]®and 
[M]* it is only required to choose proper i and j and compute 
Eqn. 3. 18 and 3.19. Moreover if all the elements have same shape 
and siae it is enough if we compute them once and use it for all 
the elements. This reduces the computation involved by an order of 


magnitude. 
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The equation (3.1) of motion can be rewritten as 


■’[ 

Oj 

- 

2 

it) 

M ‘ 


^ 1 









1 ) 


0 


The elemental matrices are assembled to get the global 
matrices. The boundary conditions are applied to global matrices 
to yield a real eigen value problem. This is solved for various 
values of C to get eigen values (frequencies) and mode shapes. 


3.5 IB<DaJ[NlDAlR\? C<DQ^IDDUD<&IHS : 

From Floquet’s theory the boundary condition for the 
considered unit cell can be written in general as (figure ^3 ) 
follows. If L,]R,B,T subscripts denote left, right, bottom and top 
boundaries of the unit cell respectively, the generalised 
displacements { V' ) along the boundary Viave the following 
relations . 


{ V = e ' t t V' 

{ V* = e ^ (3.21) 

Where k . k are the wave vectors in x, , x„ directions 

i 2 12 

respectively. 

It may be noted that that for free wave propagation in an 
infinite medium the nodal forces are zero. Thu.s for the 
equilibrium of unit cell the nodal forces are given by 

{ F^ } = - { F^> ; { F^ } = - e'^Vv^ { F^} (3.22) 
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Thus from Eq. (3.21) it can be seen that U , V , U , V_ can be 

R R T T 


expressed in terms of U^, , U^, and can be eliminated from 

the set of unknowns as follows. 




to 


[”°]J ('} 


where 


(’i 


V' 


V' 




{F}^ {0>^ -{F}^ eVa 


1 ’ 


(3. 23) 


The stiffness and mass matrices can be partitioned as follows: 


[ K 




Sym 


Kbx 


KtT 

Kbk 


(3.24) 



[ M f = 


Mlb 

Mas 


Sym 


Mli It 


(3.2&) 
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Further simplification is done using the relations Eq,3.24 and 
3.25 . The procedure is illustrated by considering the right and 
left boundary relationships. 

The right hand boundary degrees of freedom when expressed as 
a function of left hand boundary degree of freedom the system of 
equations simplify as 







c 

V + K,. K.. 


a 

¥ 

I 









r 

2 

40 




C 

V'n 





G 







(3.26) 

where )U = ik d 

Similarly the relationship between top and bottom nodal 
degrees of freedom can be used resulting in an expression similar 
to equation 3.26. Then Eqn. 3.20 can be written as, 



where: 


r 

{V'^} {w \ 

11 

1 Jr 

_ ^ 


'"i= 

and K and M matrices are real function of complex variable! le 


/j . 
^2 


Thus Equation 3,27 represents a complex eigen value problem, 
the complex quantities coming from the boundary conditions. . It can 
also seen that the matrices K and W appear as complex conjugates 
and hence are Hermit ian. 

As it is difficult to solve large eigen system with complex 
elements, it is required to convert them into system of real 
matrices. The simplification procedure is given by yjilkinson [34] 
and by J.V Sastry [35] . 

Let 


K ( - K ^ + i K ^ ; M - M i 


{ V ) = { + i (3.28) 

where r, 1 denote real and imaginary parts, Substituting Eqn.3.28 
in Eqri.3.27 the Eqn. 3.27 take,s the form 



2 


- a> 



i m’- 


y; *" + iy/ 


= 0 (3,29) 


Separating the real and imaginary parts and combining the two sets 
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of equations, Eq. 3.24 can be expressed as 



K" -iK' “ 

2 

-iM’' 



r 

r 



~ to 






k' K" 


m' If 



i 

V' 


— , 


- 








J 


L. J 


Since the matrices [ K 3 and [ M } are Hermitian , 

[ kM = - [ 3 =-[»']■' 

Thus equation 3.29 represents a real symmetric eigen value 
problem, of rank 4N where 'N ’ is the number of independent nodes 
belonging to the unit cell. The equation 3.30 can be solved for 
natural frequencies ’c*.' 'and eigen vectors { V' } by specifying the 
values of k^ , and d^ , d^ . 

3.6 CODiFiUTIEiR 0[HtPILE[MI£tNirA'irD<P[N : 

The global stiffness and mass matrices are obtained by 
assembling the elemental matrices. These elemental matrices are in 
turn obtained by calculating stiffness and mass sub— matrices for 
any pair of nodes in an element and then grouping them. 

Assembling and storage of these lai'ge matrices are effected 
using profile storage or skyline system [36], In order to apply 
the boundary conditions the boundary nodes are given the same node 
number for those of which the degrees of freedom are related. A 
special connectivity is given for the assembly purpose. Another 
connectivity and node numbering system is also provided to 
calculate the shape functions, mesh generation and calculation of 
the jacobian etc... 



3A 


Subspace iteration technique detailed in Ref36 is used for 
obtaining lowest few eigen values and eigen vectors. Because of 
the free-free boundary conditions for the cell , zero eigen values 
appear correspondiiig to the possible rigid body inodes. Hence a 
shift is applied for the eigen values [1] in order to over come 
tlie zero eigen values. 
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CHAPTER 4 

RESULTS AND DISCUSSION 

In this chapter the examples for analysis and the results 
obtained are presented. A discussion on the. results obtained is 
also presented. The following cases are solved to obtain 
dispersion relations (frequency spectrum) and mode shapes as a 
function of wave number C • 

4.1 Ei^AIHIPILIES : 

Two types of materials viz., a homogeneous isotropic material 
and composite material are investigated. V?ave propagation both in 
the plane of the plate and in the thickness plane (plane stress 
and plane strain case) is studied. An isotropic material is 
considered in order to validate the analysis procedure. Two types 
of fiber reinforced composite materials are considered for 
investigating the free wave propagation in composite materials. 
The materials of fibers and matrix are so chosen that the ratios 
of fibre and matrix stiffness are different. Properties of 
materials considered here are presented in Table 4.1. 

4.2. FD[ND¥E Efl-EIHE1N7 IMESIH : 

The unit cell size Is chosen to be ^ 

discretiaed using 8-noded isoparametric elements. The cell is 
divided into 16 elements, 4 in each direction (Fig 4.1.). The 
total number of nodes in the unit cell including those on the 
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Fig, 4,1 Finite element mesh for the unit cell 
in an infinite medium. 



Pig. 4,2 Finite element mesh for the unit cell 
in composite medium 

4,2a For propagation in plane 

4,2b For propagation in Xg - plane 
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right and top boundaries is 96. But if the nodes on the right and 
top boundaries are considered as those belonging to the adjacent 
cell but interacting with the considered cell, the number of 
nodes then are 75}belonging to the cell. 

In composite materials, the unit cells for propagation in 
the plane of fibers (plane stre.ss case ) and in the thickness 
plane (plane strain case ) are different and.aht Illustrated in 
Figs. 4,2 a.b, . From the Fig. 4. 2 a.b the volume fraction of fibres 
and matrix are as follows*. 

(a) Plane stress case ; Vj ~ 0.5 

(b) Plane strain case ; =0.3 

4,3 C.<&lN^EiR<EEtNCE CIRDFEtROOtN ; 

In the subspace Iteration the following convergence 
criterion for eigen values is used : 

x"” -x." 

V t 

where denotes the i*"*^ eigen value after the iteration 

and jp is the tolerance limit specified. 

In the present analysis two values of tolerance were tried. 
The tolerance limit of 10~® and 10“*® were provided for the eigen 



values . 
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4.4 OS<DTrtROtPD€ tHAlTE[RaAtL 

An isotropic material vjith E= 30 GPa and u -0.3 and p - 700 
kg/m^ Sixteen elements were used to discretize the unit cell. 
For this case ten lowest eigen values and corresponding eigen 
vectors are obtained for various values of wave number ranging 
from (0.001-2.00). These are obtained for both plane stress and 
plane strain cases. The results for plane stress case are 
presented in Table 4.2. and Fig 4.3. The eigen frequencies are 
normalised with respect to the third lowest frequency at = 
0.001. The normalised frequencies for plane stress case are 
given in Table 4.2 and the frequency spectrum is plotted in Fig 
4,3. 

Fi'om figure 4.3 and table 4.2 following points may be noted. 

1) The frequency spectrum is symmetric about C = 1.0 . This may 
be attributed to the fact that at C = 1 the half wavelength is 
equal to the unit cell dimension d. 

2) There are two types of branches in the frequency .spectrum 
vis., the acoustical branch and the optical branch [29]. 

In the Brillouin sense they are defined based on the way 
individual particles move in the direction of wave propagation. 
If all the points (Fig. 4. 4) within half wavelength in the 
direction of propagation are having same orientation with respect 
to direction of propagation the corresponding branches are called 
as acoustical branches. If the particle orientation is different 
at every point then it is said to belong to optical branch. 
Acoustical branches pass through the origin whereas the optical 
branches do not. 
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Table 4.^ Normalized Frequencies for an isotropic material 
^Propagation in t . *“ x plane) 



wave number zeto 


Fig. 4,3 Frequency Spectrum for an isotropic infinite 
medium for wave propagation in direction 




4 ^ 



acoustical mode. optical mode, 


Fig, 4, A Acoustical and optical modes of propagation. 



3) The acoustical branches have the same slope ( cko / dC ) 
indicating' that the group velocity is constant. This confirms the 
non dispersive effect of the isotropic medium. 

4) The optical branches in contrast, are having different group 
velocities for different C implying dispersive phenomena. 

It may be pointed out that the lowest branch corresponds to 
the shear wave propagation and the next lower one to the 
longitudinal wave in the medium, 

4.5 IRE5itia.7S-C<0MtP<DSDTE PiftTERDfttLS ; 

Two types of fiber reinforced composites are considered for 
analysis. Table 4.1 gives the material properties for the same. 
For the first one composite material the ratio E, / E = 3.276 

f- rn 

and p. J p =4 . 

f 

For the second material the ratios respectively are 20.1 and 
2.1 respectively which corresponds to Glass/Epoxy composite. Ten 
lowest eigenvalues and corresponding eigen vectors are obtained 
for various values of and varying from 0 to 2. 

Initially attempt were made to obtain the results using 
tolerance of 10’'®(Eqn. 4.1). Results showed inconsistancies and 
were not accurate. Hence tolerance has been increased to 10 
which necessitated the use of double precision aritlimetio. With 
these the results sho'wed consistency and also the convergence of 
the eigen values improved appreciably. 

Fig. 4, 5 shows rate of convergence of the eigen values with 
the eigenvalues. From the figure 4.5 it is evident that 
eigenvalues converged Quickly and after about 40 iterations there 



eigen values 





eigen values 


b 



Fig, 4.5 Convergence rate of eigen values with nmber of 
jacobi iterations, for 
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is very lit. tie defference between tbe values in subsequent 
iterations. Hence increasing the number of Iterations has little 
effect. Hence in the present analysis the .o of iterations was 
limited to 45. 

It also appeared that convergence of eigen values depends 
upon the range of elements in stiffness and mass matrices. Because 
of different material properties in different elements within the 
cell the stiffness and mass matrices tend to become ill 
conditioned. This results in increase of the number of Iterations 
and use of double precision for composite materials. 

4.5.1 COfhlP<DSD¥E [MATTlElRDAfl,- 1 ; 

Results for the composite material with E,./ E = 3.275 is 

I T17 

presented in figures 4.6a,b,c,d and tables 4. 3, 4. 4, 4. 5,4.6 and 
4.9. 

4. 5. 1.1 tPiR<0(PftC^/SgC<D[N A(L<D(NC EfHE FDtBEfRS : 

Fig 4.6a and table 4.3 gives the dispersion relations for 
the pre.sent ca.3e. Comparison of these results with those of 
isotropic material (Fig. 4,3) shows the_ following; 

1) Acoustical branches for thi.s material is almost same as to 
that of the isotropic material up to <;’^= 0.75. When C^approaches 
unity slight dispersion is evident, 

2) Optical branches indicate larger dispersion, 

3) In composites also there is symmetry in frequency spectrum 
about 1.0 
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Table 4.3 *» Normalized frequency for Compisite material < 3.275) 

Propagation along direction 




4 ■ 5. 1 . 2 [P[R<D1PAGA¥D<D1N [PE{RFE{mi>nCQJIl.ft[R irc> IFDIEIERS : 


Two cases of propagation perpendicular to fibers are 
considered ; 

a) In the plane of plate ( plane )- plane stress case 

b) In the thickness plane ( plane ) - plane strain case 

Fig. 4, 6(b) and Table 4.4 illustrate the former case while 
Fig 4.6d and Table 4.6 the second case. 

From Fig. 4. 6(b) and Table 4.4 the following points may be 
noted , 

1) Disper. 5 ive phenomenon is predominant even for comparatively 
small values of = 0.25. This is in contrast to the behavior in 
the direction of fibers ( Fig. 4. 6). Even acoustical branches show 
considerable disper.sion. 

dw 

2) The slope of the branches zero at ? =1.0 ( d = X = 

1.0) indicates the group velocity at =1.0 is zero. This can be 
attributed to the scattering of waves as the wave length approach 
the size of unit cell. 

3) Values of frequencies are appreciably less for propagation in 
direction as compared to those for direction. This can be 

attributed to the difference in stiffness of the medium in both 
the directions ( anisotropy). As the composite is more stiff in 
fiber direction the values of frequencies are more for 
propagation along the fibers. 

Fig. 4. 6(d) and Table 4.6 show the frequency spectrum for 
harmonic wave, propagation in the plane. 

Following obser'vations can be made with reference to Fig. 


4.6(d) and Table 4.6. 
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Fig, 4,61), Frequency spectrum for composite material 
(E^/E^ » 3.5L't5)for wave propagating perpendicular 
-to the fibers but in the same plane of fibers. 
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Fig.4,6'cl Freqency spectrum for wove propogotion in composite 
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1) There are some differences between the frequency spectra for 
wave propagation in plane stress and plane strain cases even if 
the direction of propagation is perpendicular to fibers. The 
difference between the first and second acoustical branches at ? 
= 1.0 is large for plane stress as compared to plane strain case. 
This is because the lowest frequency at = 1.0 for plane strain 
case is nearly twice that of the lowest frequency in plane stress 
case but the next mode value are appreciably not different. 

2) There is a additional optical branch appeax'ing In this case. 

3) In this case also the group velocity is zero at ? - 1.0 

4) There is a distinct change in the slope of the branches at 
?=0.5 

4.5.1. 3 lPtR<l>IPA<SAirDOtN [DDIRECTTOOD^ ORDEPOTID 45^-’ ITO FDtBEIRS : 

Fig. 4.7(c) and Table 4.6 show the frequency spectrum for the 
waves propagating at 45*^ ( ) to the fibers and in the 

same plane The following observations can be made. 

1) In the case of acoustical branches there is little difference 

between longitudinal and shear waves up to ? = 0.25. Then the 

difference increases with further increase in 

2) Slope of the acoustical branch.es is zero at ? -0.75 

indicating that the group velocity is zero at that point. 

3) In the optical branches the behavior is more complex. Some 

branches indicate that the group velocities at.? =0.25> ? =0.5 , 

? =0.75 and 1 are different where as some have zero group 
velocity at ? = 0 . 25 , 0.5, 0.75 
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Fig, 4,6c, Wave propagating with ^’2 

the x^- axis), E^/e^j = 20,1, 
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4) VJhen compared to frequency spectrum lor propagation along and 
perpendicular to fibers ( Fignres 4.6 (a),(b)) The absolute 
values of eigen values for lower two modes of circular waves 

have intermediate values between those for propagation along 
and perpendicular to the fiber direction. 

4.5. ^ lEO-ASS EIFOS;g'tf C(DlhlPOSOH[E: P-liElDDW ; 

A gla.ss fiber reinforced composite material was also 
investigated. Figs. 4.7 (a.b) and Tables 4.7 and 4.8 present the 
results for wave propagating in this medium. Two cases considered 
here are propagation along the fibers and 45*^ to the same, 

4.5. ^.1 lPR<PIPft<£A¥0<D{h FUE FDIBEIRS ; 

Fig. 4.7(a) and Table 4.7 gives the frequency spectrum for 
the glass fiber composites for propagation along the fibers. From 
the graph the following points can be observed. 

1) For the lowest acoustical branch for the corresponding shear 
waves there Is no dispersive effect till = 0.75 and further 
increase in ? the group velocity gradually decreases to zero at 

n = 1 . 0 . 

2) For the second acoustical branch corresponding to 
longitudinal wave there is a perceptible change in group velocity 
at f = 0.25 itself then gradually reduces to zero at = 1.0 

3) When compared to acoustical branches the optical branches 
show more complex pattern. But for all the optical branches also 
the group velocity is zero at = 1.0. 

4 ) All the optical branches show distinct change in group 
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Fig, 4,7a, Frequency spectrum for glass - epoxy composite 
medium « 20.1) For wave propagating along the 

fibers . 
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velocity at ? =0.5 

4. 5. 2.2. OF COIRCOiftlLIR MAWS; 

For propagation of circular wave ( ?^= ') the frequency 

spectrum is presented in Pig. 4.7b axid Table 4.8. From the 
propagation of waves at 45^ to the fiber orientation appears to 
be similar to that of the propagation along the fibers. 
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CHAP1*ER 5 

CONCLUSIONS 


Finite el enieni:, analysis of harmonic wave propagation studies 
in composites materials was carried out for two dimensional wave 
propagation. Concepts of unit cells and Floquet theory were used 
in modelling the behavior of the medium, he eigen value problem 
obtained by finite element formulation was solved for ten lowest 
eigen frequencies as a function of wave number C . The frequency 
spectrum (dispersion relations) were obtained for various cases of 
propagation in the plane of plate and in thickness direction with 
values of K varying from 0- 2.0 . 

Following conclusions can be drawn from the analysis carried 

out . 

1) The behavior of composite materials for harmonic wave 
propagation is quite complex when compared to that of isotropic 
materials. 

2) The ratio of stiffness of individual constituents in the 
medium is an Important factor affecting the dispersion phenomena. 

3) The dispersion phenomena in a composite medium is dependent 
on the mode of propagation and the wavelength /frequency of 
propagation. As the half wavelength approaches the dimension of 
the unit cell, dispersion is pronounced and at d = g the group 
velocity Is aero indicting scattering of waves within the medium. 

4) There is symmetry of frequency spectrum about C = 1- This may 
be attributed to the fact that fiber and matrix materials are 
symmetrically situated about the centre of the unit cell. 



6A 


£>) The diape, rsion effect depends on the direction of 
propagation. The dispersion effect is pronounced in the direction 
perpendic\.il ar to the fibers . 

6) The longitudinal and shear waves are comparatively less 
affected by the dispersive nature of the medium than the higher 
modes . 

7) Finite element method can be a suitable technique to study the 
wave propagation in composite materials. 



65 


BIBLIOGRAPHY 


1. Zienkiewicz.O.C, , " The Finite Element Method " , McGraw-Hill 
Book Company, Third Edition 1977. 


2. Osaki, T and Kimpara I , ; "Evaluation of Defects in 

Unidirectional Fiber Composite by Elastic Wave Propagation", 
Proc. of the Fourth Conf, on Composite Materials, Jime 27 - 
29, 1988. 


3. Oshima T. etal .; "Stress Wave Propagation In a Rectangula 
Composite Beam" , Proc. of the Fourth Japan -US Conf . on 
CompositeMaterials , June27-29, 1988. 


4. Navi , P. , "Waves and Vibrations in Elastic Fibre - Composite 
Medium and Plates ", Vibrations and Stress Analysis ,pp. 

135 - 143. 


5. Kohn W, J.A. Krumhansl and E.H. Lee : "Variational Methods 
for Dispersion Relations and Elastic Properties in Composite 
Materials", Trans. ASME, Jl. Applied Mech. , Vol . 39, June 

1972, pp. 327 -- 336. 


6. Dong, S.B. and Nelson R B, '• "On Natural Vibrations and 
Waves in Laminated Orthotropic Plates", Trans. ASME, Jl. 
Applied Mech. , Vol. 239, Sept. 1972, pp. 739 - 745. 


7. Yang, W H and Lee E.H. : "Modal Analysis of Floquet Waves in 
Composite Materials ", Trans. ASME, Jl. App. Mech., June 
1974, pp. 429 - 433. 


8. Richard B. Nelson and Parvia Navi ■ “Harmonic Wave 

Propagation in Composite Materials", J. of Acoust. Soc. of 
Am. , Vol. 57 , No. 4 , April 1975 , pp. 773-781. 


9. Golub, G.H., Jemming L. Yang W.H, , " Waves in Periodic 
Structured Media", J.Of Comp. Phy. , Vol 17, pp 349-357 1976 


10. Minagawa, S. et al. : "Finite Element Analysis of Harmonic 
Waves in Layered and Fiber Reinforced Composites", Int. Jl. 
for Num. Methods in Engg. , Vol, 17, 1981, pp. 1335 - 1353. 


11. Minagawa, S.et al. : "Finite Element Analysis of Harmonic 
Waves in Layered and Fiber Reinforced Composites", J. Comp. 



66 


Struc., Vol 19 (1-2), 1984, pp 119~28 

12. Dal/ta, S.K., et al. ; "Wave propagation in laminated 
composite plates" , Jl. Acoust.Soc of America. No. 63 (6) 

June 1968 


13. Mai, A.K. and Y. Barcohen : "Stress Waves in Layered 
Conipo.site Laminates", Proc. of the Fourth Japan - US 
Conference on Composite Materials, June, pp. 27 - 29, 1986. 


14. ToshiukiOshima , et al . , " Stress Waves in a Rectangular 
Composite Beam " , Proc., Foi'th Japan U.S., conf Comp. Mat., 
Jan 27-29 1986 


15. Masanory Koshiba " Finite Element Analy,sis of Wave Scattring 
in an Elastic Plate Wave Gaide ", IEEE Trans Sonic and 
Ultasonicsa Vol Su 31 No, 1 Jan 1984 


16. Joseph L. Rose : "Ultrasonic Wave Propagation Principles in 
Composite Material Inspection ", Materials Evaluation / 43 / 
April 1985, pp. 481 - 483. 


17. Rose, J.L., A.S.D. Wang and E.W. Leska : "Wave Profile 
Analysis in a Unidirectional Graphite Epoxy Plate", J. of 
Composite Materials, Vol. 8, Oct., 1974, pp. 419 - 425. 


18. Tauchert, T.R. and Gazel.su ; "An Experimental Study of 
Dispersion of Stress Wave in a Fiber Reinforced Composite", 
Trans. ASME, J. of App. Mech. , Vol. 79, Mar., 1972, pp. 95 - 
102 . 


19. Munson, D.E. and K.W. Schuler : "Steady Wave Analysis of 
Wave Propagation in Laminar and Mechanical Mixtures", J. of 
Composite Materials, Vol. 5, July 1971, pp. 286. 


20. Hemann and George Y. Baaklini "The Effect of Stress on 
Ultrasonic Pulses in Fiber Reinforced Composites", SAMPE J., 
July / Aug. 1986, pp. 9. 


21, Teti, R and G. Caparino : "NDE of THick GFRP Composites 
through Ultrasonic Wave Form Detection" , Developments in the 
Science and Technology of Composite Materials, ECCM -* 3, 3rd 
European Conference on Composite Materials, 20 - 23 .March, 
1989, Bourdeaux, France, pp. 793 - 800. 


22. Kinra, V.K. and V. Dayal : "A New Technique for Ultrasonic - 



67 


Non destructive Evaluation of Thin Bpecimens" , Experimental 
Mechanics, Vol . 28, No. .3, Sept. 1988, pp. 288 ~ 297. 


23. Navadunsky, J.J. and Lucas, J.J., ; “ Early Fatigue Damage 
Location In Composite Materials”, J. Comp. Mat. Vol 9 Oct 
1975 pp 394-408 


24. Rose, J.L., and Carsonm. , : " Ultasonic Nondestructive Test 
Procedure For The Early Detection of Fatigue Damage and The 
Prediction of Remaing Life " , Mat . Eval. , Apr 1980 pp 
27-34 


25. Jones, R.M. '• "Mechanics of Composite Materials”, 1975. 


26. Pei , Chi Chen Wang A.S.D. , " COntrol Volume Analysis of 
Eleastic Wave Front in Composite Materials”, J.Comp. Mat 
Vol 4 Oct 1970 pp 444-461 


27. Tewary, V.K. : "Mechanics of Fiber Composites", Wiley 
Eastern Ltd., 1978. 


26. Lawrence J .Broutman and Richard H.Krock Composite 

Materials " , Vol. 2, " Mechanics of Composite 

Materials”Edited by Sendeokyj , G.P., Academic Press, NY, 
1974. 


29. Kittlel.,C., “Intoduction To Solid State Physics ” Willey 
Eastern Edition, Newyork, 3 Edition 1968 . 


30. Kolsky, H. : "Stress Waves in Solids", Dover Publishers Lnc. 
NY, 1963. 

31. Humar, J.L. : "Dynamics of Structures", Prentice Hall Ltd., 
1969. 

32. Chueng. ,Y.K. , and Yeo, ,M. P. , " A Practical Introduction To 
Finite Element Analysis", Pitman Publishing Ltd, 1978 


33. Cook, R.D. • "Concepts and Applications of Finite Element 
Analysis", John Wiley and Sons., 1981. 

34. Wilkinson, " The Algebriac Eigen Value Problem , Oxford: 
The Clerendon Press 1976 



68 


3&. Sastry, J.V. ’’In PLane Vibration of Layered Rings on 
Periodic Radial Supports Using Finite Element Method 
M.Tech Thesis Dept. Mech Engg. , I.3.T., Kar2pur,13B6 


3G . Bathe, K.J. and E.L. Wilson : ’’Numerical Methods in Finite 
Element Analysis” , PHI Pvt. Ltd., 1967. 


37. Reddy, J.N. ; ” An Introduction to the Finite Element 
Analysis ”, 1986. 



AtoH'iez 


7h 


0<D*/97 

TiIit book ii to be returned on the 


date last stamped. 




vJav 



